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^ ' Abstract. In this paper we study the shape differentiability properties 

, of a class of boundary integral operators and of potentials with weakly 

' singular pseudo-homogeneous kernels acting between classical Sobolev 

I spaces, with respect to smooth deformations of the boundary. We prove 

^N) , that the boundary integral operators are infinitely differentiable without 

loss of regularity. The potential operators are infinitely shape differen- 
, tiable away from the boundary, whereas their derivatives lose regularity 

near the boundary. We study the shape differentiability of surface dif- 
ferential operators. The shape differentiability properties of the usual 
strongly singular or hypersingular boundary integral operators of inter- 
est in acoustic, elastodynamic or electromagnetic potential theory can 
' then be established by expressing them in terms of integral operators 

$_( . with weakly singular kernels and of surface differential operators. 
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fundamental solution, surface differential operators, shape derivatives, 
Sobolev spaces. 



1. Introduction 

Optimal shape design problems and inverse problems involving the scattering 
of time-harmonic waves are of practical interest in many important fields of 
applied physics including radar and sonar applications, structural design, bio- 
medical imaging and non destructive testing. We develop new analytic tools 
that can be used in algorithms for the numerical solution of such problems. 

Shape derivatives are a classical tool in shape optimization and are 
also widely used in inverse obstacle scattering. In shape optimization, where 
extrema of cost functions have to be determined, the analysis of iterative 
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methods requires the study of the derivative of the solution of a scatter- 
ing problem with respect to the shape of the boundary of the obstacle. An 
explicit form of the shape derivatives is required in view of their implemen- 
tation in iterative algorithms such as gradient methods or Newton's method 
[H m [24]. By the method of boundary integral equations, the shape analysis 
of the solution of the scattering problem with respect to deformations of the 
obstacle is obtained from the Gateaux differentiability analysis of boundary 
integral operators and potentials with weakly singular, strongly singular, or 
hypersingular kernels. An expression of the shape derivatives of the solution 
can then be computed by taking the derivative of its integral representation. 
This technique was introduced for the Dirichlet and Neumann problems in 
acoustic scattering by Potthast dU [52] and applied to the Dirichlet problem 
in elastic scattering by Charalambopoulos [T] in the framework of Holder con- 
tinuous and differentiable function spaces. More recently these results were 
exploited in acoustic inverse obstacle scattering to develop novel methods in 
which a system of nonlinear integral equations has to be solved by a regular- 
ized iterative method [HI [T31 [H] . 

An extension of the technique to elasticity and electromagnetism re- 
quires the shape differentiability analysis of the relevant boundary integral 
operators. More generally, we are concerned in this paper with the Gateaux 
differentiability of boundary integral operators with strongly and weakly sin- 
gular pseudo-homogeneous kernels acting between classical Sobolev spaces, 
with respect to smooth deformations of the boundary considered as a hyper- 
surface of M.'^ with d S N, d > 2. This family of integral operators covers 
the case of the single and double layer integral operators from the acoustic 
and the elastic scattering potential theory. The differentiability properties of 
the hypersingular boundary integral operators can then be established by ex- 
pressing them as products of integral operators with weakly singular kernels 
and of surface differential operators. In return, however, we have to study the 
shape differentiability of surface differential operators. The electromagnetic 
case presents a specific difficulty: The associated boundary integral operators 
act as bounded operators on the space of tangential vector fields of mixed 
regularity TH~ ^ (Jivp^ r). The very definition of the shape derivative of an 
operator defined on this energy space poses non-trivial problems. This is the 
subject of the second part of this paper [3] where we propose an analysis 
based on the Helmholtz decomposition [3] of TH~ = (divr, F). 

This work contains results from the thesis [17] where this analysis has 
been used to construct and to implement shape optimization algorithms for 
dielectric lenses, aimed at obtaining a prescribed radiation pattern. 

The paper is organized as follows: 

In Section [5] we describe the family of pseudo-differential boundary in- 
tegral operators and potentials that we consider. We use a subclass of the 
class of pscudo- homogeneous kernels introduced by Nedelec in his book [20] . 
Main results on the regularity of these operators are set out. In Section[31 we 
define the notion of shape derivative and discuss its connection to Gateaux 
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derivatives. We also recall elementary results about differentiability in Frechet 
spaces, following ideas of [SUS] and notations of pS] , 

Section H] is dedicated to the shape differentiability analysis of the inte- 
gral operators. We discuss different definitions of derivatives with respect to 
deformations of the boundary and compare them to the notions of material 
derivatives and shape derivatives that are common in continuum mechanics, 
see Remark 14.11 We prove that shape derivatives of the boundary integral 
operators are operators of the same class, that the boundary integral opera- 
tors are infinitely shape differentiable without loss of regularity, and that the 
potentials are infinitely shape differentiable away from the boundary of the 
obstacle, whereas their derivatives lose regularity in the neighborhood of the 
boundary. A main tool is the proof that the shape differentiability of the inte- 
gral operators can be reduced to the one of their kernels. We also give higher 
order Gateaux derivatives of coefficient functions such as the Jacobian of the 
change of variables associated with the deformation, or the components of 
the unit normal vector. These results are new and allow us to obtain explicit 
forms of higher order derivatives of the integral operators. A utilization for 
the implementation of higher order iterative methods is conceivable. 

The shape differentiability properties of usual surface differential op- 
erators is given in the last section. Again we prove their infinite Gateaux 
differentiability and give an explicit expression of their derivatives. These 
are then applied to obtain the derivatives of hypersingular boundary integral 
operators from acoustic, elastic and electromagnetic potential theory. 

Notice that our shape differentiability analysis is realized without re- 
striction to particular classes of deformations of the boundary, such as it is 
frequently done in the calculus of variations, namely restriction to deforma- 
tions normal to the surface as suggested by the structure theorems for shape 
derivatives [SI [9l [21] , or consideration of radial deformations of star-shaped 
surfaces [2l[l3l[l2]. 



2. Pseudo-homogeneous kernels 

Let O denote a bounded domain in M.'^ with d > 2 and let ff^ denote the 
exterior domain R"^ \ fi. In this paper, we will assume that the boundary F 
of ri is a smooth closed hypersurface. Let n denote the outer unit normal 
vector on F. 

For a domain G C M.''' we denote by H''{G) the usual L^-based Sobolev 
space of order s S M, and by H^^^{G) the space of functions whose restrictions 
to any bounded subdomain B oi G belong to ^{B). 

For any t G K we denote by H*(r) the standard Sobolev space on the 
boundary F. The dual of i?*(F) with respect to the scalar product is 
iJ^*(F). Vector functions and spaces of vector functions will be denoted by 
boldface letters. 
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Q\a\ 

For a = (tti , . . . , Ud) £ N'^ and 2 = (zi , . . . , z^) G M"* we denote by — — 
the linear partial differential operator defined by 

Q\a\ Qai gad 



dz°' dz^' dz^" 

where \a\ = ai + ■ ■ ■ + ad- For m G N, the total differential of order m, a 
symmetric m-linear form on R'^, is denoted by D™. 

The integral operators we consider can be written in the form 

]Cru{x)= k{y,x-y)u{y)ds{y), x G F, (2.1) 



where the integral is assumed to exist in the sense of a Cauchy principal value 
and the kernel k is regular with respect to the variable ?/ G F and pseudo- 
homogeneous with respect to the variable z = x — y G M''. We recall the 
regularity properties of these operators on the Sobolev spaces -ff*(F) for all 
t G R, available also for their adjoint operators 

IC*riu)ix) = k{x, y - x)u{y)ds{y), x G F. (2.2) 

We use a variant of the class of weakly singular kernels introduced by Nedelec 
in m pp. 168ff]. More details can be found in [3 [H IH H HSl US • 

Definition 2.1. The kernel G{z) G '^°°(M'' \ {0}) is said to be homogeneous 
of class —m for an integer m > if 

(i) for any a G there is a constant Cq, such that for all z G R'* \ {0} 



we have 



5; 



-G(z) 



Q\a\ 

(ii) for any a G N'' with lal — m, the function — — G(z) is homogeneous 

dz" 

of degree — (d — I) with respect to the variable z, 

(iii) D™ G(z) is an odd function of z. 

Remark 2.2. (i) The number — m in this definition is not the order of homo- 
geneity of the kernel, but related to the order of the corresponding pseudo- 
differential operator defined on the d — 1-dimensional manifold F. 
(ii) Our condition (iii) is stronger than the vanishing condition in Nedelec's 
original definition, but it is easier to verify, and it is satisfied for the classical 
integral operators we will be considering. 

Definition 2.3. The kernel k{y, z) defined on F x {R'^ \ {0}) is said to be 
pseudo-homogeneous of class —m for an integer m such that m > 0, if the 
kernel k admits the following asymptotic expansion when z tends to 0: 

k{y,z)= brn+j{y)GU,{z), (2.3) 

j>0,i 
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where for j — 0, 1, ... the sum over £ is finite, &^_|_, belongs to "^""(r) and 



rn+j 

^m+j is homogeneous of class —{m + j). 

In (|2.3p . one can also consider coefficient functions of the form brn+j{x, y) 
with X = y + but using Taylor expansion of such coefficients at z = 0, we 
see that this would define the same class of kernels as with (12.31). 



Example 2.4. (Acoustic kernels) Let k S C \ {0} with Im(K) > and d — 2 
or d — 5. The fundamental solution 

^-H^^\k\z\) whend = 2 

— -— when d = 3 

47r|z| 

of the Helmholtz equation Am + k^u = in R'' is pseudo-homogeneous of 
class —1. Its normal derivative ^^^(y) priori pseudo-homogeneous 

of class but one can show that in the case of smooth boundaries it is a 
pseudo- homogeneous kernel of class -1. 
Indeed one can write 

The first term is homogeneous of class —1, the second term is smooth and 
for j > 3 the j-th term is homogeneous of class —(1 -I- j). The double layer 
kernel has the expansion 

g^fa{n,z) = n(,).V^G.(.,z) = (n(y).z) ^ ^ ^ 

One can prove that the function g{x,y) = n[y) ■ [x — y) behaves as |x — 
when z = X — y ^ Q (see for instance [IHl P- 173]). We refer to example 14.111 
for a proof using a local coordinate system. 

Example 2.5. (Elastodynamic kernels) Let a; G K and d = 2 or d = 3. Denote 
by p, /i and A the density and Lame's constants. The symmetric fundamental 
solution of the Navier equation — /xAu — (/z -I- A)Vdivii — pw^u = 0, given 

by 

Ge{Ks,Kp, z) = - [GaiKs: z) ■ l^d + HcSS {Ga{K.s,z) - GaiUp, z)) 



with Kg = ^\Jji ^-nd Kp = \+2n ' pseudo-homogeneous of class —1 . 
The traction operator is defined by 

du 

Tu — 2u— h A(div u)n + an A curlu. 

on ^ ' 

The double layer kernel (TyGei^s, i^p,x ~ y))~^ is pseudo-homogeneous of 
class 0. The index y of Ty means that the differentiation is with respect 
to the variable y. Notice that TyGe{Ks, Kp, x — y) is the tensor obtained by 
applying the traction operator Ty to each column of Ge{Ks, Kp,x — y). 
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For the proof of the foUowing theorem we refer to [3S] . 

Theorem 2.6. Let k be a pseudo-homogeneous kernel of class —m. The asso- 
ciated boundary integral operator JCr given by (|2.ip is linear and continuous 
from -ff*(r) to _ff*"'"™(r) for all < G M. The same result is true for the adjoint 
operator /Cf . 

The foUowing theorem is estabhshed in [7]. 

Theorem 2.7. Let s G M. Let k be a pseudo-homogeneous kernel of class — m. 
The potential operator V defined by 

V{u){x) = j^k{y,x~y)u{y)ds{y), x £ M"* \ F (2.4) 

_i_ _i_ 1 S I TTh I 

is linear and continuous from H^{T) to 2 u Hi^^ ^ (^'^)- 

3. Some remarks on shape derivatives 

We want to study the dependence of operators defined by integrals over the 
boundary F on the geometry of F. This dependence is highly nonlinear. The 
usual tools of differential calculus require the framework of topological vector 
spaces which are locally convex at least, a framework that is not immediately 
present in the case of shape functionals. The standard approach consists in 
representing the variations of the domain VI by elements of a function space. 
We consider variations generated by transformations of the form 

X ^ X r{x) 

of point X in the space IR'', where r is a smooth vector function defined in 
the neighborhood of F. This transformation deforms the domain VI into a 
domain Vlr with boundary F^. The functions r are assumed to be sufficiently 
small elements of the Frechet space X — '^°°(F,M'^) in order that (1 + ^") is a 
diffeomorphism from F to 

F^ = (I + r)F = [xr ^x + r{x)\x G F} . 

For e small enough we set 

S°°(0,e) = {r e ^°"{tX): doo(0,r) < e] , 

where c?oo is the distance induced by the family of non-decreasing norms 
(II • ||fe)feeN defined by 

||r|U = sup sup |D™r(a;)| . 

0<m<fc a:eR<i 

Consider a mapping F defined on the set {F^; r G B°°{Q, e)} of bound- 
aries. We introduce a new mapping 



B°^{0,e) 3r^ Tr{r) = F(F^). 
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We define the shape derivative of the mapping F through the transformation 

Tb x + i{x) e R'^ by 

dF[r; Hm = li„, MtQ - MO) 

t^O t t^Q t 

if the hmit exists and is finite. The shape derivatives of F are related to the 
Gateaux derivatives of Jt (see 0121]). 

Fix ro G i?°°(0,£). Following the same procedure, one can construct 
another mapping -^r,.,, defined on the family of boundaries 

{(I + r')(r.J; r'eB-(0,£')}- 

Notice that J"r.„(0) = F{Tr,) = J"r(ro) and J"r.„((r - ro) o (I + ro)-i) = 
F{Tr) = Mr). 

3.1. Differentiability in Frechet spaces: elementary results 

Frechet spaces are locally convex, metrisable and complete topological vector 
spaces on which the differential calculus available on Banach spaces can be 
extended. We recall some of the results. We refer to Schwartz's book [33] for 
more details. 

Let X and y be Frechet spaces and let U he & subset of X. 

Definition 3.1. (Gateaux semi-derivatives) The mapping / : J7 — > 3^ is said 
to have a Gateaux scmiderivativc at tq G J7 in the direction of ^ G A" if the 
following limit exists in y 

t^o t at\t=o 

Definition 3.2. (Gateaux differentiability) The mapping f : U ^ y is said to 
be Gateaux differentiable at tq G C/ if it has Gateaux semiderivatives in all 
directions ^ G A" and if the mapping 

X3^^df[ro;^]ey 

is linear and continuous. 

We say that / is continuously (or '^^-) Gateaux differentiable if it is 
Gateaux differentiable at all tq G U and the mapping 

U xX 3df ■.iro,0^df[ro;C](^y 

is continuous. 

Remark 3.3. In the calculus of shape derivatives, we usually consider the 
Gateaux derivative at r = only. This is due to the result: If Jr is Gateaux 
differentiable on B°°{0,e), then for all ^ (£ X we have 

dFr[ro;^] = dF[Tr„;^ o (I + ra)-^] = dTr^^, [0; f o (I + ro)-^]. 
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Definition 3.4. (Higher order derivatives) Let m G N. We say that / is 
(m + l)-tinies continuously (or Gateaux differentiable if it is 

Gateaux differentiable and 

U5r^d^f[r;^i,...,U] 

is continuously Gateaux differentiable for all m-tuples {£,1 , ■ ■ ■ , £,m) S X™. 
Then for all tq G C/ the mapping 

A-^+i 9 (Ci, . . . , U+i) ^ Ci, • ■ • , U+i] e y 

is (m + l)-linear, symmetric and continuous. We say that / is "^""-Gateaux 
differentiable if it is "^'"-Gateaux differentiable for all to G N. 

Proposition 3.5. Let f : U ^ y be "t^"^ -Gateaux differentiable. Let us fix 
roeU and^e X. We set -f{t) = f{ro + t£). 

i) The function of a real variable 7 is of class ^™ in the neighborhood 
of zero and 



l^'^Ht) = d^\^J('o + tO = d"^f[ro;t_^]- (3.2) 

m times 



ii) We use the notation 

d 



We then have 



1 0?n 
p—1 l<'ii<---<2p<m 



(3.3) 



Thus the knowledge of — — / Wo ; ^] suffices to determine the expression 
ofd^/iro;^,...,^- 

Proposition 3.6. Let f : U -^ybe -Gateaux differentiable. Let us fix 
ro G U and £ E X with £ sufficiently small. Then we have the following 
Taylor expansion with integral remainder : 

fc=i •^'J 

The chain and product rules are still available for ^"-Gateaux differ- 
entiable maps between Frechet spaces. 
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4. Shape differentiability of boundary integral operators 

Let Xr denote an element of Tr and let be the outer unit normal vector 
to r^. When r = we write no — 'H- We denote by ds{xr) the area element 
on Tr- 

In this section we want to establish the differentiability properties with 
respect to r G i?°°(0,e) of boundary integral operators /Cr^ defined for a 
function Ur G H^(Tr) by: 

(ICr.Ur) (Xr) ^ J kr{yr,Xr - yr)Uriyr)ds{yr), Xr (4.1) 

and of potential operators Vr defined by: 

(VrUr) (x) = J fcr(yr, X - yr)Ur{yr)ds{yr) , X £ fi,. U fij;, (4.2) 

where K € "^"^ (r^ x {R"^ \ {0})) is a pseudo-homogeneous kernel of class 
— m with m 6 N. 

We point out that we have to analyze mappings of the form r i~->- J-r (r) 
where the domain of definition of J-rif) varies with r. This is the main dif- 
ficulty encountered in the calculus of shape variations. We propose different 
strategies according to the definition of the mapping J^r- 

(i) A first idea, quite classical (see HH [H]), is that instead of studying 
mappings r i— )■ J-r{r) where ^r(^) = Ur is a function defined on the boundary 
Tr, we consider the mapping 

r 1-^ Ur o (1 + r). 

Typical examples of such functions Ur are the normal vector rir on Tr and 
the kernel kr of a boundary integral operator ICr^ (see Examples l2.4l and l2.5p . 

To formalize this, we define the transformation ("pullback") which 
maps a function Ur defined on Tr to the function o (I + r) defined on T. 
For all r G B°°(0, e), the transformation is linear and continuous from the 
function spaces '^''(Tr) and H*'{Tr) to '^'^(F) and i?*(F), respectively, and 
admits an inverse. We have 

{TrUr){x) — Ur{x + r[x)) and {T^'^u){Xr) = u{x) . 

(ii) Next, for linear bounded operators between function spaces on the bound- 
ary, we use conjugation with the pullback t^: Instead of studying the mapping 

B°°(0,£) 3r^^ J-r(r) =/Cr, G ^ (iJ^(F,), /i'^+'"(F,)) 

we consider the mapping 

B^{0,e)3r^TrJCr,T-^ G ^ (iJ"(F), iJ^+™(F)) . 

We have for u G iJ'*(F) and x G F: 

{TrlCr,T-^){u){x) =y kr{y+r{y),x+r{x)-y~r{y)) u{y) Jr{y)ds{y), (4.3) 
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where Jr is the Jacobian (the determinant of the Jacobian matrix) of the 
change of variables on the surface, mapping x to x + r{x) G F^. 

(iii) The third case concerns potential operators acting from the boundary to 
the domain: 

Each domain $7 is a countable union of compact subsets: 51 = IJ Kp. For 

pSN 

ah p e N, there exists Sp > such that Kp C H VLr- Thus, instead of 

rGS(0,ep) 

studying the mapping 

B°°{Q,e) 3 r ^ Frir) = Vr ^ ^ (H'{Tr),H'+"'+^{nr)) 
we can consider the mapping 

B°°{Q,ep) 3r^ P,r-i G ^ (i?^(F), (ifp)) . 

We have for u G 77* (F) 

{Vr{r)T;r^){u){x) = J kr{y + r{y),x - y - r{y)) u{y) Jr{y)ds{y), x G Kp. 

(4.4) 

Then passing to the hmit p — J> cxd we can deduce the differentiabihty proper- 
tics of the potentials on the whole domain il. We use the analogous technique 
for the exterior domain fl'^. 

In the framework of boundary integral equations, these approaches were 
introduced by Potthast [HJ (55] in order to study the shape differentiability 
of solutions of acoustic boundary value problems. 

Remark 4.1. In continuum mechanics, when the deformation x > r{x) — 
ro{x) + tS,{x) is interpreted as a flow with initial velocity field S,{x), one 
frequently considers two different derivatives of functions Ur defined on O^. 
The material derivative iir is computed by pulling Ur back to the reference 
domain fJ, thus by differentiating r H> r,.Ur = Uro{l-\-r). The shape derivative 
u'^{x) at a point x is defined by differentiating Ur{x) directly. At r = the 
difference between the two derivatives is a convection term: 

"0 "o + C • Vuo • (4.5) 
This is easily seen from the definition of the material derivative 

Ur{x) = d{Tu)[Q]S]{x) = ut(,{x + ti{x)) = d u[Q; ^]{x) + £,{x) ■ V Uq{x) . 

dt t=o 

Relation ()4.5|) can be used to compute the shape derivative from the simpler 
material derivative, see jTJ] for an application. 

In this terminology, the derivatives of boundary functions and operators 
in (i) and (ii) above would be analogous to material derivatives, whereas the 
derivatives of potentials in (iii) correspond to shape derivatives. Instead of 
formally defining the terms "material derivative" and "shape derivative" , we 
prefer here to explain in each instance precisely which Gateaux derivative 
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is meant. We want to emphasize, however, that the shape derivatives of so- 
lutions of electromagnetic transmission problems can be obtained by using 
the three kinds of derivatives defined above. This will be explained in detail 
in Part II of this work. The construction is based on an integral representa- 
tion of the solution of the transmission problem by potentials, the densities 
of which are solutions of boundary integral equations with operators of the 
type studied here. Thus the mapping from the given right hand side to the so- 
lution is a composition of boundary integral operators, inverses of boundary 
integral operators, and potential operators. By the chain rule, its derivative 
is then obtained by composing boundary integral operators, their inverses, 
and potential operators with derivatives of type (i), (ii), and (iii) above. The 
same structure gives the shape gradient of shape functionals that are defined 
from the solution of the transmission problem. In this case, also adjoints of 
the boundary integral operators have to be differentiated. This poses no new 
problem, because adjoints of operators with quasi-homogeneous kernels have 
quasi-homogeneous kernels, too. 



4.1. Gateaux difFerentiability of coefficient functions 

For the analysis of the integral operators defined by (|4.3p and (|4.4p . we first 
have to analyze coefficient functions such as the Jacobian of the change of 
variables T 3 x i-^ x -\- r{x) S F^, or the normal vector on F^. 

We use the standard surface differential operators as described in detail 
in [20]. For a vector function v G '^'=(K'^, C'') with k eW,we denote by [Vv] 
the matrix the i-th column of which is the gradient of the i-th component of 
V, and we write [Dv] = [Vv]^ . The tangential gradient of a scalar function 
u e 'r'=(F,C) is defined by 

Vru V{t|r - (Vu|r • n) n, (4.6) 

where u is an extension of u to the whole space R'^. For a vector function 
u G '^^{T, C^), we again denote by [Vrw] the matrix the i-th column of which 
is the tangential gradient of the i-th component of u and we set [Dp u] ~ 
[Vr«]^. 

We define the surface divergence of a vector function u G '^'''(F, C^) by 

divr u — divM|r — ([VM|r]w • n) — divu|r — ■ ' (4-7) 

where u is an extension of u to the whole space M'*. These definitions do not 
depend on the choice of the extension. 

The surface Jacobian Jr is given by the formula Jr — Jacr(I+?') = \\wr\\ 

with 

Wr = cof(l + Dr|r)n = det(I + D r|r)(I D r|r)"^^n, 

where cof (A) means the matrix of cofactors of the matrix A, and the normal 
vector rir is given by 
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The first derivative at r = of tliese functions are well known, we refer for 
instance to Henrot-Pierre [S] . Here we present a method that allows to obtain 
higher order derivatives. 

Lemma 4.2. The functional J mapping r S B°°(0,e) to the Jacobian Jr S 
'^'^{TjM.) is -Gateaux differentiable and its first derivative at tq is given 
fori e •=^'^{T,W^) by 

dj[ra,(\ = Jro(T"ro divr^Jr-^O)- 

Proof. We just have to prove the '^°°-Gateaux differentiability of 

W : B°°{0,e) 3 r ^ Wr ^ cof(I + Dr|r)n e ^°°(r). 

We use a local coordinate system. Assume that F is parametrized by an atlas 
{Oi, 4>i)i<i<p then Tr can be parametrized by the atlas (O^, (I + r) o 4>i)i<i<p. 
For any x S F, let us denote by ei(x), 62(0;), ... , ed-i{x) a vector basis of the 
tangent plane to F at a;. A basis of the tangent plane to F^ at a; + r{x) is 
then given by 

ei(r, a;) = [(I + T>r){x)\ei{x) for i — 1, . . . , d — 1. 

Notice that for i = 1, . . . , d — 1 the mapping B°°{0,e) 3 r H> ei{r) e 
^°°(F,M'') is '^""-Gateaux differentiable. Its first derivative is de^ro;^] = 
[D^]ei(ro), and higher order derivatives vanish. We have 



d-1 

A ej(r,a;) 



where the wedge means the exterior product. Since the mappings r n- ei{r), 
for i = 1, . . . , d — 1 are "^""-Gateaux differentiable, by composition the map- 
ping W is, too. We compute now the derivatives using formulas 
Let f G ^°°(F,R'^) and t small enough. We have at tq € B°°(0,£) 



d-1 

A (I 



Dro + tDi)e^{x) 



d-1 

A e^ix) 



To simplify this expression one notes that 

[D^(a:)]e,(x) = [B ^ixMl + B ro){x)]-'[{l + B ro){x)]e,{x) 

= [De(x)][D(I + ro)-\x + ro{x)mi + Bro){x)]e,{x) 
= KD(Tr^^O(a;)]e»(»^o,a;) = [tto Dr^o (T"ro^O(2;)]ej(ro, a;). 

Now given a {d x d) matrix A we have 

d-1 d-1 

^ • • • A Ci-i X Aei A Cj+i A • • • = (Trace(^)I - A'^) /\ e^. 
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Thus we have with A = [r,.„ Dp,.,, (t^^^O] and Bq = I, Bi{A) = Trace(^)I-AT 

Jro ( yrn 

[i?i(A)e]>V(ro), 
[B,„(A)e]W(ro) 



>V(ro) 

aw 



i=l [ITT' ^ ^j- 

for 1 < m < d — 1 



9™W 
It follows that 



[ro, ^] = for all m > d 
1 dW 



dr ^""^'^^ ||W(ro)|| dr 



ro,e]-W(^o) 



9W 
i9r 



[ro, C] • TroWro = Jro (Tro divr,„ ('^ro^O) • 



□ 



From (#) we deduce easily the Gateaux differentiability of r i-^ r^rir- 

Lemma 4.3. The mapping Af from r G B°°{0,e) to Trtir ~ rir o {I + r) S 

'^°°(r,M'') is "1^°° -Gdteaux-differentiable and its first derivative at tq is de- 
fined for ^ e 'r°°(r,R'*) by: 

|^[ro,?] = -KVr.„(r-/e)]AA(ro). 

Proof. Using the preceding proof, we find 



1 



/-I 



W(ro)||3 V 5^ 



□ 



To obtain higher order shape derivatives of these mappings one can use 
the equahties (#) and 

WTrllrW = 1, 



dr" 



for all m > 1. 



For example, we have at r = in the direction ^ e "^^"^(F, M.'^): 

d 7 dJ\f 

^[0,^] = divr^ and ^[0,^] = -[VrC]". 
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Using Proposition [Xni we obtain 



0,fi,6] = -Trace([Vr6][Vr6])+divra-divr6 + ([Vr6]"- [Vr6]«). 
e that Trace([Vr6][Vrei]) = Trace([Vr6][Vr6])- 
^ 0,^1,6] = [Vr6][Vrei]"+ [Vra][Vr6]n- ([Vrei]"- [Vr6]n)n. 



In the last section we give a second method to obtain higher order derivatives 
using the Gateaux derivatives of the surface differential operators. 

Remark 4.4. The computation of the derivatives does not require more than 
the first derivative of the deformations ^. As a consequence for hypersurfaces 
of class it suffices to consider deformations of class '^'^+^ to conserve 

the regularity of the Jacobian and of the normal vector by differentiation. 

4.2. Gateaux differentiability of pseudo-homogeneous kernels 

The following theorem establishes sufficient conditions for the Gateaux dif- 
ferentiability of the boundary integral operators described above. 

Theorem 4.5. Let p e N. We set (r x T)* = {(x, y) e T x T; x ^ y}. Assume 
that the following two conditions are satisfied: 

1) For all fixed (x, y) G (F x F)* the function 

f : B°°{0,s) C 

r 1^ kr{y + r{y),x + r{x)-y-r{y))Jr{y) 

is '^P'^^ -Gateaux differentiahle. 

2) The functions {y, x — y) ^ fifa)iy, x ~ y) and 

{y,x- y) ^ d''f[ro,£,i, . . .,^i]{y,x- y) 

are pseudo-homogeneous of class —m for all ro G B°°(0,e), for all I — 
l,...,p+l and for all ^i, . . . , ^p+i € 'r°°(r, R''). 
Then for any s G M the mapping 

B°°{0,e) ^(i/''(F),i/''+'"(r)) 
is 'i^P -Gateaux differentiahle and 

(F (r,./Cr,T7^) [ro,Ci, • ■ ■,ip\u{x) = j d^ f[rQ.Xi, ■ ■ .,^p]iy,x- y)u{y)ds{y). 

Proof. We use the linearity of the integral and Taylor expansion with integral 
remainder. We do the proof for p = 1 only. Let tq G B°° (0,£), ^ G '^°°(F,E'') 
and t small enough such that tq + G B°°{0, e). We have 

df 

f{ro + t^,x,y) - f{ro,y,x- y) = t—[ro,(,]{y,x- y) 



+ t^ 



/ V - A) [ro + {y, X - y)dX. 



We have to verify that each term in this equality is a kernel of an operator 
mapping iJ'*(F) to iJ'*+™(F). The two first terms in the left hand side are 
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df 

pseudo- homogeneous kernels of class — m and by hypothesis -^[ro,^] is also 
a kernel of class — m. It remains to prove that the operator with kernel 



I (l-A)0[ro + Ate,e](x,y)dA 



acts from jy*(r) to i/''+'"(r) with norm bounded uniformly in t. Since 
92/ 

-tt^Itq + At^,^] is pseudo- homogeneous of class —m for all A e [0,1], it 

suffices to use Lebesgue's theorem in order to invert the integration with 
respect to the variable A and the integration with respect to y on F. 



{^-^)g^[ro + Xt^,iKx,y)dX ] u{y)ds{y) 



(1-A) 



H»+'"(r) 



< 



suPag[o,i] 



[ro + Xt£„£,]{x,y)u{y)ds{y) 



< C\\u\ 



(r) 



We then have 
1 



t 



f{ro +t^,x,y)u{y) ds{y) 
df 



dr 



f{ro,x,y)u{y) ds{y) 
ro,£,]{x,y)u{y)ds{y) 



+ t 



S_]{x,y)dX u{y) ds{y). 



We pass to the operator norm limit t ^ and we obtain the first Gateaux 
derivative. For higher order derivatives it suffices to write the proof with 
(i^/[ro, ^1, . . . , ^fc] instead of /. The linearity, the symmetry and the continuity 
of the first derivative are deduced from the corresponding properties of the 
derivatives of the kernel. □ 

Now we will consider some particular classes of pseudo-homogeneous 
kernels. 

CoroUjiry 4.6. Assume that the kernels are of the form 

kriVn Xr ~ yr) — G{Xr ~ yr) 

where G G '^°"(R'^\{0}) is a pseudo-homogeneous kernel of class —m, m E N, 
which does not depend on r. Then the mapping 

B°°(0,£) ^ J^(i/*(F),i?*+™(F)) 
r TrlCr^T^^ 

is "^"^ -Gateaux differentiable and the kernel of the first derivative at r = is 
defined for <^°°(F,R'^) by 

df[0, = i^x) - C(y)) • VG(x -y) + G{x - y) divr C(y). 



16 



Martin Costabcl and Frederique Le Louer 



Proof. For fixed {x,y) G (F x F)*, consider the mapping 

/ : B°°(0, e) 3 r ^ /(r, x, y) - G{x + r{x) -y- r{y))Jr{y) £ C. 



By Theorem 133] we have to prove that r ^ f{r) is "^""-Gateaux difFerentiable 
and that each derivative defines a pseudo-homogeneous kernel of class — m. 
> Step 1: 

First we prove that for fixed {x,y) G (F x F)* the mapping r H> f(r,x,y) 
is infinitely Gateaux differentiable on B°°{0,e). By Lemma the mapping 
r I— Jr{y) is infinitely Gateaux differentiable on B°°{Q,e), the mapping r i-> 
x + r{x) is also infinitely Gateaux differentiable on i3°°(0, e) and the kernel G 
is of class on IR''\{0}. Being composed of infinitely Gateaux differentiable 
maps, the mapping r i— f(r,x,y) is, too. 
l> Step 2: 

We then prove that each derivative defines a pseudo-homogeneous kernel of 
class — m, that is to say that for all p € N and for any p-tuple (^i, . . . , ^p) the 
function 

(x,y) dPf[ro,^i, . . .,^p]{x,y) 

is pseudo-homogeneous of class — m. By formula p.3p . it remains to write the 
gp 

proof for the function — — fWo,^] with ^ G '^°°(F,R'*). The Leibniz formula 
orP 

gives 



QP P /p\ Ql 

f[r,M^,y)=Y.[i)Q^l 
n \ 



QP-I 7 

{G{x + r{x) - y ~ r(y))}[ro,C]^^[ro,e](2/)- 



Since — rb'O:^] G "^""(F,!!), we have to prove that 

fjj^p^i 

Ql 

{x, y)"^ {^(^ + "^i^) f{y))} [ro, £,] 

defines a pseudo-homogeneous kernel of class — m. We have 

Ql 

-g^{G{x + r{x) - y - r{y))}[ro;^] 

= D' G[x + roix) - y - My); ^x) - ^(y), • . • , - ^y)]. 

By definition, G{z) admits the following asymptotic expansion when z tends 
to zero: 

N-l 

G{z) = G^{z) + G,n+ji^) + G,n+N{z) (4.8) 

where Gm+j is homogeneous of class — (m+j) for j = 0, . . . , A^— 1 and Gm+N 
is of arbitrary regularity. Using Taylor expansion, the following result is easy 
to see: 

Lemma 4.7. Let the kernel Gm{z) be homogeneous of class — m and ^ G 
^°°(r,R'^). Then the function 

[x,y-x) ^ D'Gmix + ro{x) - y - ro{y);£,{x) ~ S,{y), . . . ,£.{x) - ^{y)] 
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is pseudo-homogeneous of class —m. 

By taking derivatives in the expansion ()4.8p we conclude that 



TT-r {G(x + r(x) — y ~ ^{y))} [tq; is pseudo-homogeneous of class —m too. 
or 

This ends the proof of the coroUary. □ 

Theorem 4.8. Let s G M. Let G{z) be a pseudo-homogeneous kernel of class 
— {m-\- 1) with m G N. Let us fix a compact subdomain Kp ofQ. Assume that 
for all r G B°°{0, Bp), we have kr{yr, x — yr) — G{x — yr). Then the mapping 

^ (^H''-i{T),H''+"^{Kp)'j 
r ^ VrT-^ 
is infinitely Gateaux differentiable and 

dP{Prr-')[ro,(i,...,^p]u{x) 

dP {G{x -y- r{y))My)} [ro,a, • • ■ ,^pHy)ds{y). 

r 

Its first derivative at r = in the direction ^ G '^°°(r,IR'') is the integral 
operator denoted by "P^^-* with kernel 

-ay) ■ y'G{x ~y) + G{x - y) divr i{y). 

The operator V'^^^ can be extended to a continuous linear operator from 
H'-^T) to H'+"'{n) and H'i+J^(Jf^). 

Proof. The kernel and its higher order derivatives are of class "^"^ on Kp. 
Writing f2 as an increasing union of compact subsets, we can define a shape 
derivative on the whole domain $7. Let us look at the first derivative: The 
term G{x ~ j/)divrf(2/) has the same regularity as G{x — y) when x — y 
tends to zero wheareas f (y) • VG(x — y) loses one order of regularity. As a 
consequence, since the kernel is of class — (m -|- 1), its first derivative acts 
from H"-^ (r) to and iJf +" (fF) . □ 

Remark 4.9. We conclude that the boundary integral operators are smooth 
with respect to the domain whereas the potential operators lose one order of 
regularity at each derivation. We point out that we do not need more than 
the first derivative of the deformations ^ to compute the Gateaux derivatives 
of any order of these integral operators. 

Example 4.10. (Acoustic single layer potential) Let d = 2 or d = 3 and s G M. 

We denote by the single layer potential defined for Ur G H^iTr) with the 
fundamental solution Ga of the Helmholtz equation (see Example 12. 4|) 



5'^Ur(a;) = / Ga{K,x - yr)uriyr)ds{yr), X eR \ 
Let its trace on 

V^Urix) = I Gain,X- yr)Ur{yr)ds{yr) , X G F^. 



18 



Martin Costabcl and Frederique Le Louer 



Since Ga is pseudo- homogeneous of class —1, the mapping 

B°°(o,e) ^ ^(i?"(r),ij"+i(r)) 

is infinitely Gateaux differentiable. The mapping 

r ^ Tr^^T^i 

is infinitely differentiable and its first derivative at r = can be extended to 
a Unear continuous operator from H^{T) to H'^^^{Q) U (^'^)- 

Similar results can be deduced for the elastic single layer potential. 

Example 4.11. (Acoustic double layer kernel) Let d ^ 2 or d ~ 3 and t G M. 

We denote by 13^ the boundary integral operator defined for Ur G i/*(r,.) by 

D^Urix) = J TlriXr) ' VGa(K, J/r " )Wr (j/r )ds (j/r ) ■ 

The mapping 

S°°(0,e) ^(i7*(r),iJ*+i(r)) 

is ^°°-Gateaux differentiable . 
Indeed the mapping 

B°°(0,e) 9 r 5(r, a;, y) = (T^n^)(a:) • (a; + r(a;) - y - r{y)) 

is "^"^ Gateaux differentiable and by using a local coordinate system (see 
[22] ) we prove (when d = 3) that the Gateaux derivatives behaves as \x — y\'^ 
when X — y — 0. We use the same notations as in the proof of Lemma Fix 
a; G r and set gxir,y) = g{r,x,y). We have that g^ G 'rf'="(B°°{0,e) x r,M). 
If r is parametrised by the atlas {Oi, 4'i)i<i<p then when a; G = 0^(0^) nF 
we can write x = 0i(?7i,?72) where (^yi,??!) G ^i- The tangent plane to F at 
a; is generated by the vectors ei(x) = f^('?f,f?f) and 62(0;) = ^{Vi,V2)- 
Thus gx{r, (f>i{r]i,'i]2)) has the expression 

(I + D r ) f| (r;f , 77f ) A (I + Dr) |^ jr^f , T^f ) 
I (I + Dr) f| (,7f , ryf ) A (1 + Dr) f| (,7? , r;f ) | ' 

• {{l + r) o(f,i{ijf,r]^) - (l + r) 00,(771,772)) 
Using Taylor expansion we have when y — x 

9x{r,y) =0 + J:>g^{r)[x;y - x] + g^{r)[x;y - x,y - x] + . . . 

Writing gx{r) = {gx{r) o t/ij ) o t/ir ^ , we have for all r G B°°{0, e) that 
Tlga^ir) = D(^^^,,3)(5^(r) o0,) oD^ri. 

By straigthforward computations we obtain that D(^j^^ jj^){gx{r) o i^j) = for 
all r \2l\ . Thus by differentiation with respect to r we prove that g^ {r, y) and 
all its Gateaux derivatives behaves as \x — ?/|^ when a; — y — > 0. 
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5. Shape differentiability of surface differential operators, 
application to hypersingular boundary integral operators 

Many classical hypersingular boundary integral operators can be expressed 
as compositions of boundary integral operators with pseudo-homogeneous 
weakly singular kernels and of surface differential operators. Such repre- 
sentations are often used in the numerical implementation of hypersingular 
boundary integral operators. Here we use these representations to study the 
shape derivatives of hypersingular boundary integral operators. To this end, 
in addition to the shape derivatives of the weakly singular integral boundary 
integral operators as studied in Sectional we need to determine the Gateaux 
derivatives with respect to deformations of the surface differential operators 
acting between Sobolev spaces: The tangential gradient is linear and contin- 
uous from iJ*+i(r) to H*(r), the surface divergence is linear and continuous 
from H*+i(r) to i?*(r). 

Example 5.1. (Acoustic hypersingular kernel) Let k £ C with Im(K) > and 
d = 3. The hypersingular kernel is the normal derivative of the double layer 
kernel. We have 

dn{x) dniyf''^'''''"^^ 

= -n{x) ■ n{y)AGa{K,x - y) + n{x) ■ curF (V^Gq(k,x - y) A n{y)). 

When d — 2, for a scalar function (p the term —\7(p A n is the arc- length 

derivative Using integration by parts with respect to the variable y and 
ds 

that for a scalar function v and a vector a G M'' it holds n ■ curl(i;a) = 
— (Vu An) ■ a we obtain for a scalar density u 

[ n{x)-curr {vy{Ga{K,x-y)) An{y))u{y)ds{y) 

r 

= {\7-{Ga{K,x-y))Anix)) ■ (uiy)) A n{y)) ds{y) . 
Finally we have 

- J^{V'fGa{K,x-y) An{x)) ■ {\/ru{y) An{y))ds{y). 

A similar technique can be applied to the elastic hypersingular boundary 
integral operator using integration by part and Giinter's tangential deriva- 
tives (see [TTIITB]). 

Lemma 5.2. Let d — 3 and T be a closed orientable surface in M"^. The 
tangential Giinter derivative denoted by Ai is defined for a vector function 
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V e 'ri(r,c3) by 

d 

M.V — iVv — (div v) ■ lR3)r2. = -tt—v — (div v)n + n A curl v. 

on 

(i) We set n = {nk)i<k<^ and My = {'mjk)i<j,k<j,- We have 

nijk = nk{y)^ ^jiy)-^ = -mkj. 

oyj dyk 

(ii) For any scalar functions u,u in "^^(F, C) and vector functions v,v in 
"^^(F, C^) there holds the Stokes formula 



(mjku) ■ uds = — J u ■ (rrijku) ds and J (Mv) ■ v ds — + J v ■ (Mv) ds. 

(5.1) 

Example 5.3. (Elastic hypersingular kernel) Let w G M and d — 3. Denote 
by p, /i and A the density and Lame's constants. The hypersingular kernel is 
defined by 

H{x,y) ^ T^{TyGe{K,x - 

where Ge is the fundamental solution of the Navier equation and T is the 
traction operator defined in Example 12.51 First of all we rewrite the operator 
Tu as 

Tu = 2pMu + (A + 2^)(divM)n - fj,n A curlw. (5.2) 

Then we apply the operator Ty in the form (|5.2p to the tensor Ge{ns, 'tp, x—y). 
It follows 

{TyGe{Ks,Kp,x - y)y 

= 2fi{MyGe{K, X - y))^ - {n{y) A curlj^ Ga{Ks,x - ?/)Ir3)^ 

+ ^^'^ (niy) ■ divy Ge{Ks,Kp, x - y))^ . 

dlYy Ge{Ks,Kp,X - y) 

= {VyGaiKs,X - ?/))''' + ^VyAy(Ga(Ks,X - y) - Ga{Kp,X - y)) 



-^{VyGaiKp,x - y)y 

n{y)Ac\iv\yGa{ns,x-y)lf^3 = ^TWy - + ri{y) • divy^ Ga{K.s,x-y)l^3 

In virtue of the property (i) in Lemma 15.21 we can write 
[n{y) AcMr\yGa{ns-,x-y)ljfj^ = [-My - q^^^^^ Ga{Ks,x- y)ljg,3 

+ {n{y)-VlGaiK„x-y)y 
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Collecting the equalities we obtain 

{TyGe(Ks,Kp,x - y))~^ 

= 2n{MyGe{K,x -y))^ + (^ gyf(^) +My^ Gains, x - y)lR3 

+ Vy{Ga{np,x - y) ~ Ga{Ks,x- y)j ■ n{yY . 

By integration by part and using the properties (i) and (m) of Lenima l5.2l we 
obtain that 

(TyGeiKs, Kp, X - y))~'^u{y) ds{y) = 2/i J Ge{Ks, Kp, x-y)Myu{y) ds{y) 

f d 

Ga{Ks,x — y)A4yu{y) ds{y) + J Q^^y^ Gai^s, x — y)u{y)ds{y) 



V 



r 



y{Ga{Kp,x- y) - Gains, X- y)j (n(y) • uiy)) dsiy). 



The kernel of the last term in the right hand side is pseudo-homogeneous of 
class —2. Thus applied to this term yields a pseudo-homogeneous kernel of 
class —1. Similarly to {TyGeing, Kp, x — y)) , the kernel T^Geins, Kp, x — y) 
can be rewritten in terms of products of weakly singular kernels and the 
Giinter derivative Mx- Now we apply the operator to the kernels of the 
second and third terms on the right hand side in the form 

T^u = (A + /i)n( divj; u) + ^ g^^^^ + u. 

We obtain 

^^^Ga(Ks,X-2/)-lK3j+(A + /i)n(.T)-Vj^^^Ga(Ks,X-2/) 



- Tx\Gains,x - y) ■ 1^3 1 = Gains, x - y) ■ Ir3 

- ^j.MxiGaiKs,x - y) - Irs) - (A + /i)n(a;) ■Vx^GaiKs,x - y) 
We use the equality 



d_ 
' dn{y) ' 



^^"Er77T'^''('^«'^ ~ y) = MyVxGaiKs,x - y) - niy)AyGains,x - y) 
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and Lemma [5.21 to show that (see [TT] pp. 52) 
f d 

I n{x)-Vj GaiKs,x ~ y)u{y)ds{y) 

Jr dn{y) 

- / n{x) -yjGains.x ~ y)Myu{y)ds{y) 

Finally we have 

T^{TyGe{Ks, Kp,x- y)yu{y) ds{y) 

= 2fi J [T^Ge{Ks,Kp,x - y)]Myu{y)ds{y) 



^1 9n(x)9n(^) ^°("-^'^)"^^^^^(^^ 



r 

f d 

/ o t . Ga{Ks,x-y)Myu{y)ds{y) 
-^lMx J GaiKs,x-y)MyU{y) ds{y)+^Mx J g^^^^ GaiKs,x-y)u{y)ds{y) 
+ J TxVy{Ga{K,p,x - y) - Ga{Ks,x - y)^{n{y) ■ u{y)) ds{y) 
+ „2(A + ,W.,/G.,».,.-.,(nM. 



r 

We see that the boundary integral operator with either the acoustic 
hypersingular kernel or the elastic hypersingular kernel are operators of order 
+ 1 on the Sobolev spaces H*{T) for i e M. Using the integral representations 
above, the differentiability properties of these operators can be deduced from 
the knowledge of the differentiability properties of the surface differential 
operators. Following the same pullback procedure as in Section|H the analysis 
of the hypersingular integral operators is finally reduced to the analysis of 
the mappings 

r ^ Tr divr,, t^^ . 

Indeed, the Giinter derivative can be rewritten in terms of these two differ- 
ential operators: 

M.V = (Vrf — (divr v) ■ iRs^n. 
The results are established in the following theorems. 

Theorem 5.4. The mapping 

g : B°°{0,e) J^(i7*+i(r),H*(r)) 
r TvVr^ 
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is '^^'^ -Gateaux differentiable and its first derivative at ro is defined for ^ £ 
'^°°(r,]R'') by 

dg[ro,^]u = -[g{ro)^]g{ro)u+{g{ro)u-[g{ro)m{ro))Af{ro). 

Remark 5.5. Note that we can write dAf[ro,^] = —[g{ro)^]N'{ro). Since the 
first derivatives of Af and g are expressed in terms of TV and g, we can obtain 
the Gateaux derivatives of ah orders recursively. 



Proof. In accordance with the Definition (|4.6p and Lemma 14731 to prove the 
"^""-Gateaux differentiability of g we have to prove the "^""-Gateaux differ- 
entiability of the mapping 

f : B^{0,e) 3 r ^ S^u ^ Tr (s/^uj ^ | £ ^(if*+i(r), H*(r)). 
For X G r, we have 

Tr(s7Tr^uj^ (x) = V (u o (I + r)"^) l^ {x + r{x)) 

= (I + Dr)f^;^(x + r(x)) oVS|,(x), 

and 

(I + D r)p^; {x + r{x)) = [(I + D r),, (x)] . 

The mapping g : B°°{0,e) 9 r n> (I + Dr)|j, G ^°°{T) is continuous, and 
"^""-Gateaux differentiable. Its first derivative is dgp,^] = [D^]|j, and its 
higher order derivatives vanish. One can easily see that the mapping h : r £ 
Bf ^ {x^ [g(r)]-i(x)} G <^°°(r) is also "^"^ Gateaux-differentiable and 
that we have at tq and in the direction ^: 

dh[ro,£] = -/i(ro) o dg[ra,S] o /i(ro) -h{ra) o [D^]|j, o h{ra). 

and 

d"/i[ro,ei,...,U = i-irY. (I+D'^o)"'o[T,„Dr-/e.(i)]°.-.°KDT,-ie.(„)] 

where is the permutation group of {1, . . . , n}. Finally we obtain the 
Gateaux differentiability of / and we have 

df[ro,^]u = -[f{roW{ro)u. 

Notice that this result can also be justified by using commutators : for ex- 
ample at r = in the direction , we have 

d 

where ttt = f • V. 

To obtain the expression of the first derivative of g we have to differen- 
tiate the following expression: 

g(r)u — (TrVr,,T^^u) = TVV (^Tr^^uj — (^TrTlr ■ ^T^V (''''^^"))) ''V^r 

= /(r)u-(/(r)u-A/-(r))AA(r). 
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By Lemnia l4.3l and the chain and product rules we have 

+ {f{ro)u ■ [g{ro)miro))^(ro) + {f{ro)u-M{ro)) [G(ro)e]A/'(ro) 
Combining the first two terms in the right hand side, we get 

dg[roA] = -[a(ro)^]/(ro)u+ (/(ro)u-A/-(ro)) [e(ro)W^o) 
+ (/(ro)u-[a(ro)e]A/'(ro))A/'(ro) 
= ^[g{ro)i]g{ro)u + {f{ro)u-[g{ro)(]M{ro))M{ro). 
To conclude, it sufhces to note that 

{f{ro)u ■ [g{r,)m{ro)) = (a(ro)u • [g{r,KW{ro)) . 

□ 

Theorem 5.6. The mapping 

V: B°°(0,e) ^ ^(H*+\r), i7*(r)) 
r ^ Tr divr^ TjT"'^ 

is -Gateaux differentiable and its first derivative at vq is defined for f G 
'^°°(r,]R'') by 

dV[ro,^]u = -Tmce{[g{romg{ro)u]) + {[g{ro)u]M{ro) ■ [a(ro)?K(ro)) . 

Proof. For u G H*^^(r) we have 'D{r)u = Trace([C/(r)tt]). Then we use the 
differentiation rules. □ 

Remark 5.7. (i) Since the first derivative of V is composed of g and TV and 
the first derivative of J is composed of J and P, we can obtain an expression 
of higher order derivatives of the Jacobian recursively. 

(ii) Denoting by Alp the tangential Giinter derivative on F^, the formulas 
(#) in section m can be rewritten as 

W(ro) ^,(T^„n^J, 
^[-o,C] = -J.„(-.o-^r,,(r-^e)), 

-^^[fQ,'?] = OforaUTO>(i. 

Remark 5.8. (Electromagnetic hypersinguljir kernel) Let k G C with Im(K) > 
and d = 3. The electromagnetic hypersingular operator is defined for a 
tangential density] G TH*(r) by 

C^\{x)^-- / n{x) A{c\XYr cnvr {Ga{K.,x-y)i{y)))ds{y). 
Jr 

Using the identity curl curl = — A + V div we have 



CJ(a;) = -n(x) A / (^uGaiK, x - y) ■ I 



+ -W^GaiK, x~y) divr ]j{y) ds{y). 
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This is the operator of the electric field integral equation in electroniagnetism. 
The operator is a priori an operator of order +1 on the space of tangential 
vector functions TH'(r), but it is well known that this operator is a bounded 
Frcdholm operator on the space of tangential vector fields of mixed regular- 
ity TH~ 2 (clivp, r), the set of tangential vector fields whose components are 
in the Sobolev space H~2(T) and whose surface divergence is in H~^(r). 
Therefore it is desirable to study the shape differentiability of this opera- 
tor defined on the shape dependent space TH~ = (divp, F). For this, the tools 
presented above are not directly applicable. It is the purpose of the second 
part [3j of our paper to present an alternative strategy using the Hclmholtz 
decomposition of the space TH~ ^ ((;[ivp^ r). 
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